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Active Control of Wake/Blade-Row Interaction Noise 
Through the Use of Blade Surface Actuators 

Kenneth A. Kousen and Joseph M. Verdon 
United Technologies Research Center 
East Hartford, Connecticut 06108-1049 


Summary 


This report describes a combined analytical/computational approach for controlling of 
the noise generated by wake/ blade- row interaction through the use of anti-sound actuators on 
the blade surfaces. A representative two-dimensional section of a fan stage, composed of an 
upstream fan rotor and a downstream fan exit guide vane (FEGV), is examined. An existing 
model for the wakes generated by the rotor is analyzed to provide realistic magnitudes for the 
vortical excitations imposed at the inlet to the FEGV. The acoustic response of the FEGV 
is determined at multiples of the blade passing frequency (BPF) by using the linearized 
unsteady flow analysis, LINFLO. Acoustic field contours are presented at each multiple of 
BPF illustrating the generated acoustic response disturbances. Anti-sound is then provided 
by placing oscillating control surfaces, whose lengths and locations are specified arbitrarily, 
on the blades. An analysis is then conducted to determine the complex amplitudes required 
for the control surface motions to best reduce the noise. It is demonstrated that if the number 
of acoustic response modes to be controlled is equal to the number of available independent 
control surfaces, complete noise cancellation can be achieved. A weighted least squares 
minimization procedure for the control equations is given for cases in which the number of 
acoustic modes exceeds the number of available control surfaces. The effectiveness of the 
control is measured by the magnitude of a propagating acoustic response vector, which is 
related to the circumferentially averaged sound pressure level (SPL), and is minimized by a 
standard least-squares minimization procedure. 
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1. Introduction 


Aerodynamically-induced engine noise is becoming an increasingly important subject. 
The noise produced by an engine arises from a variety of sources, but of particular concern 
are the discrete frequency tones generated through fan rotor/stator interaction. These tones 
are quite annoying to a listener, and represent a significant environmental concern [1]. The 
discrete tones appear in an engine noise spectrum as spikes at multiples of the blade passing 
frequency (BPF), whose magnitudes are significantly higher than the background broadband 
noise levels. It has been observed that noise reduction methods aimed at lowering broadband 
noise can make the perceived noise problem worse, by further emphasizing the discrete 
tones. Consequently, the turbomachinery designer is faced with the responsibility of reducing 
discrete tone noise directly. 

A major source of engine discrete tone noise is that generated within a fan stage by the 
interaction of the wakes from the fan rotor with the blades of the fan exit guide vane (FEGV). 
This interaction creates acoustic response waves, which under certain conditions propagate 
unattenuated in amplitude as they move upstream and downstream from the FEGV. The 
propagating acoustic response modes therefore provide targets for control analyses aimed at 
reducing discrete tone noise. 

In this report, we examine the noise generated by wake/blade- row interaction in a fan 
stage, and investigate one possible strategy, namely the generation of anti-sound by the os- 
cillations of blade surface actuators, for reducing this noise. Control is provided by adjusting 
the amplitudes of the control surface, or piston, motions on the FEGV blades, to produce 
acoustic response modes that are 180 degrees out of phase with the undesired response. The 
superposition of the acoustic responses due to external excitations, i.e., those associated 
with the rotor wakes, with those due to the control surface motions, results in a reduced net 
acoustic response. 

Active control of noise is an active area of research [2, 3, 4, 5, 6, 7], but applications to the 
turbofan engine noise problem have only recently begun to appear [8, 9]. The experimental 
studies in [8] and [9] suggest that the introduction of antisound as a means of active control 
of turbomachinery noise may be a promising approach. The present report provides a com- 
putational simulation of turbomachinery noise under realistic operating conditions, and its 
subsequent control. 

Section 2.1 of this report presents the analyses used to determine the vortical and entropic 
excitations that represent the rotor wakes as seen by the stator. The linearized unsteady 
analyis, LINFLO [10, 11], which is reviewed in Section 2.2, provides a means for determin- 
ing the acoustic response of the FEGV to the vortical and entropic excitations. Since the 
acoustic response of the FEGV is primarily attributable to vortical excitations, the entropic 
excitations are not employed in the subsequent control analysis. Section 3 contains con- 
trol algorithms developed to reduce the undesired acoustic responses. It is observed that if 
the number of independent control surface actuators is equal to the number of propagating 
acoustic response modes, then oscillation amplitudes of the actuators can be found that 
result in a complete elimination of the noise. Analyses are also provided, based on weighted 
least-squares minimization procedures, for reducing the noise as much as possible in the 
more general case, where the number of propagating response modes exceeds the number of 
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control surfaces. Results of the active control analyses, both for complete cancellation and 
least squares minimization, are shown in Section 4. These analyses are intended to provide 
a demonstration of the potential of active noise control through the use of blade surface 
actuators, as a method for reducing discrete tone turbomachinery noise. 
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2. Physical Problem 


A two-dimensional representation of a typical fan stage is shown in Fig. 2.1. In each row, 
upstream flow quantities are indicated by a subscript -oo, while downstream quantises are 
shown by a subscript +00. The wheel velocity of the fan rotor cascade is Vw. Total (V) and 
mean (V) velocity components are given in the figure with and without tildes, respectively. 
Where confusion between reference frames might occur, primes are used to denote quantities 
associated with the fan rotor. 

Axial and circumferential distances are represented via the Cartesian coordinate systems 
(£', r]') and (£, Tj) rigidly attached to the fan and FEGV rows, respectively. For the FEGV row, 
the mean or steady-state positions of the blade chord lines coincide with the line segments 
rj = £tan0 + mG, 0 < ( < cos0,m = 0, ±1, ±2, ... , where m is a blade number index, 
0 is the cascade stagger angle, and G is the gap vector which is directed along the 77-axis 
with magnitude equal to the blade spacing. Fan cascade quantities are defined analogously. 

The velocity triangle relationships representing the transformation of velocities from 
rotor-fixed to stator-fixed coordinates are examined in the general discussion of the wake 
analysis below, which determines the vortical and entropic excitations resulting from the 
fan wakes. Subsequently, the linearized unsteady flow analysis LINFLO is reviewed, which 
shows how the vortical excitations are used in a numerical calculation of the resultant acous- 
tic response of the FEGV. 

2.1 Wake Excitation Model 

The LINFLO analysis determines the acoustic response of a blade row to individual 
Fourier modes of a vortical or entropic disturbance imposed at the inlet to the blade row. The 
goal of the wake analysis, therefore, is to determine these Fourier modes through an analysis 
of the fan wakes and a transformation from rotor-fixed to stator-fixed coordinates. The flow 
within each wake is assumed to be parallel and aligned with the mean exit flow direction of 
the rotor, and is characterized by a given velocity profile, centerline velocity defect VT n /V e ', 
and wake half-width 8 , which are defined below. The wakes have a circumferential period 
equal to the rotor gap G' . In addition, it is assumed that across each wake, the static 
pressure P' and the total enthalpy H' T are constant. The assumption of constant static 
pressure across the fan exit implies that the fan wakes will give rise to velocity disturbances 
that are purely rotational or vortical in nature. Since the orthogonal components of the 
velocity disturbances associated with vorticity are not independent, it is sufficient to find 
the individual Fourier components of v • ejv, i.e., the component of the disturbance velocity 
normal to the inlet flow direction of the FEGV, at each multiple of BPF. These quantities 
will be found in terms of the wake shape parameters, the stator inlet flow properties, the 
rotor wheel speed, and the gap to chord ratios for each row. Finally, the assumption of 
constant total enthalpy provides a link between the disturbance velocity components arising 
from the wakes, and the associated disturbance entropy components. 

The velocity triangle for the region between the rotor and stator is shown in Fig. 2.1. 
Both the total velocity and the mean velocity from the rotor row, designated V' +oc and V' +00 
respectively, are seen to be oriented at the exit flow angle fl^oo of the rotor. The rotor moves 
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Figure 2.1: Velocity triangles for fan rotor/fan exit guide vane configuration. 

at a wheel speed of V\y in the positive rf direction. By adding the wheel velocity vectorially 
to the mean and total exit flow velocities of the rotor, the mean and total inlet flow velocities, 
V-oo and V-oo, respectively, to the stator can be found. The difference between the mean 
and total inlet velocities to the stator defines the perturbation velocity v, i.e., 

* = v;^- v;^. (2.1) 

Consequently, the perturbation velocity can be found by an analysis of excitation row quan- 
tities alone. 

In keeping with the notation of Ref. [12], the symbols and u* will be used to represent 
the components of v in directions normal and tangential to the inlet mean flow velocity of 
the stator (Fig. 2.2). From the figure, the component is found to be 

v + = V • e„ = (V^ - VU sin(n' +00 - f^) , (2.2) 

where e N is a unit vector normal to the stator inlet flow direction. 

In the rotor frame of reference, the wake velocity can be expressed as a Fourier series in 
the circumferential direction, whose zeroth component can be identified as the mean velocity. 
Consequently, v + can also be expressed as a Fourier series, whose individual components are 
given by 

OO 

v + ( t )')= X expCm/c^') , (2.3) 
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Figure 2.2: Magnification of the circled region from the velocity triangle. 

where k v = 27 r/G' is the circumferential wavenumber of the vortical excitation, G' is the gap 
to chord ratio for the rotor, and 

w n = 77, F* v-e^exp(-m«^)^', ( 2 - 4 ) 

G J K' ( 

where is a reference location in the rotor frame. 

Note that the rotor-fixed and stator- fixed coordinates are related by x' = x — Vwf, where 
x' = (CtV ) i s a position vector in the moving frame attached to the rotor, and x = (£,7/) 
is a position vector in the stationary reference frame. Therefore, the v+ components defined 
by Eq.(2.3) in the rotor frame can be expressed as 

OO 

v + = ^2 t>+ exp(mK„»;)exp(— in/c„VW*) (2-5) 

n= — oo 

in the stator frame, which demonstrates that while the disturbance velocity is steady in 
the rotating frame, it provides a set of unsteady disturbances in the stationary frame with 
frequencies u n = —tik^Vw- The interblade phase angles a n of these disturbances are found 
from the circumferential wavenumber and the gap to chord ratio of the stator row, i.e., 
<r n = nK„G. 

Wake parameters 

To evaluate the v+ components, it is thus necessary to have an expression for the rotor 
exit velocity, V' . This velocity will be evaluated by using the semi-empirical rotor wake 
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model of Majjigi and Gliebe [13], which defines the wake velocity at streamwise locations 
downstream of the rotor in terms of a shape function, a centerline velocity defect, and a 
wake half-width. The behavior of these quantities will be described below, as will their 
implementation into the analysis for the calculation of the vortical and entropic excitations 
to the FEGV. 

The rotor wakes are assumed to be aligned with the exit flow direction of the rotor, 
and to be identical from blade to blade. In terms of a Cartesian coordinate system ( T', N ') 
tangential and normal to the rotor exit flow angle, with origin at the rotor blade trailing 
edge, the velocity in a reference wake can be expressed as 

nr, N') e T , = V„'(T') (l - e r ., (2.6) 

where V^'(T') is the velocity at the edge of a wake, Vj rf (7V') is the wake deficit, defined as 
the difference between the edge velocity and the velocity in the wake, at each TV', is the 
minimum streamwise velocity across the wake, and ey is a unit vector tangential to the rotor 
exit flow direction. The quantity F d ' ef / V^ n = /(TV') is a shape function that is dependent only 
on the normal coordinate, TV'. The wakes are assumed here to be symmetric, so the minimum 
velocity occurs at the wake centerline and the quantity V^ /VJ will be hereafter referred 
to as the centerline velocity defect. Two different shape functions are usually considered, 
a hyperbolic secant function or a Gaussian function. A given wake velocity distribution is 
made unique by specifying the shape function, the centerline velocity defect, and the wake 
half-width, which is defined as the distance between the points where the wake deficit is half 
its maximum value. The shape functions are therefore given by 

/(TV') = sech ^(cosh -1 2)TV'j (2.7) 

for the hyperbolic secant profile, or 



for the Gaussian profile. The hyperbolic secant profile is used in the results presented below. 
The centerline defect VT n /V^ and the half- width 6 are given by the empirically determined 
functions [13] 


YLn = c i/4 ( 0-3675 T' + 1.95 \ 
V' D \ 7.65 T' + 1.0 ) 


(2.9) 


and 


0.2375Cp /8 T' + 0.034125 
0.357 Co /8 T'+ 1.0 


( 2 . 10 ) 


where Cp is the section drag coefficient. In these equations, the streamwise distance T' 
and the wake half-width S have been nondimensionalized by the rotor chord. A typical 
wake velocity profile, based on the hyperbolic secant shape function, is shown in Fig. 2.3. 
Fig. 2.4 shows the variation of the centerline velocity defect and the half-width with increasing 
streamwise distance from the fan rotor under consideration. Note that the decay of the 
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Figure 2.3: Typical wake profile V'(N) = V e ' - where V' c is the edge velocity, 6 is the 

wake half-width, and V^ ef {N') is the velocity deficit at each N'. 


centerline velocity defect occurs primarily in the near wake. 

After determining the velocity distribution across the wake, the Fourier decomposition 
in Eq. (2.4) must be performed on the component of this velocity normal to the inlet flow 
direction for the stator, at an appropriate axial location downstream of the rotor. The 
results of the Fourier decomposition then yield the values in Eq. (2.3). Note that care 
must be taken to ensure that all velocities are nondimensionalized with respect to the mean 
inlet flow velocity to the stator. In terms of a given centerline velocity defect, the rotor wake 
edge velocity can be shown to be 



( 2 . 11 ) 


where the mean exit flow velocity of the rotor can be expressed in terms of the stator row 
inlet flow angle and the wheel speed by 


= \fl + V|y - 2V v v sin 


( 2 . 12 ) 


Values for the tijj^s have been determined at several axial stations within the region 
between the fan rotor and FEGV. The wake profiles were determined by evaluating the 
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Figure 2.4: Variation in centerline velocity defect, V^ lin /V r e / , and wake half-width, 6, with stream- 
wise distance along a blade wake for a representative fan rotor. 

parameters 8 and n /V^ at each axial station through the empirical relations in Eqs. (2.9) 
and (2.10), after a transformation from (£',17') to ( T',N ') coordinates. Typical results for 
the variation in the magnitudes |u+|, n = 1,2, ..., 6 , of the complex amplitudes with 
increasing axial distance are shown in Fig. 2.5. Note that there is an initial decay in the 
magnitude of each component followed by a much more gradual decrease with axial distance. 
Also, the decay is more pronounced for the higher order components. 

The inviscid analysis of the FEGV imposes the vortical exitations at the inlet of the 
FEGV row. Since vortical and entropic disturbances are convected through the response 
row, their magnitudes are unaltered. Figure 2.5 demonstrates that if the axial spacing 
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Figure 2.5: Variation of the magnitudes of the complex Fourier components |u+| with axial 
distance downstream of the trailing edge of the fan rotor. 


between the blade rows is sufficiently large, the magnitudes of the Fourier components will 
not vary substantially. It is necessary, however, to choose an axial location at which the 
t>+ components will be evaluated. Common practice suggests that the leading edge of the 
response blades is a reasonable choice, and will be adopted herein. 

Entropic disturbances 


It has been assumed that the fluid pressure and total enthalpy are constant across each 
wake. These restrictions can be used to relate the entropic disturbances associated with the 
wakes to the velocity disturbances associated with the wake vorticity. The total enthalpy 
downstream of the fan rotor is given by 


m 


1 ~ 


T,+oo 


= Koo + o V +°o • VU = 7(7 - ir'PUKoo + «V' +0C . V' +00 , 


(2.13) 


where H' +00 is the static enthalpy. After linearizing this relation, setting Hj,+oo = #r,+oo> ^+00 = 
P+oo, and introducing the thermodynamic relation s' +co = 7 - 1 p+ 00 / P +00 — P+ 00 / P + we find 
that 


= -7-'(T - = -(7 - l)v;„ • (2-W) 

where A 2 = 7 Pjp is the speed of sound propagation in the mean flow. From Fig. 2.1, it can 
be seen that • v ' +00 = In addition, s' + <x> = s_oo and = u-oo- Therefore, 
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the disturbance entropy that acts as the excitation at the inlet of the FEGV is given by 

•S-co = (1 — 7)-^+oo^-°o/-^+oo (2.15) 

where CLoo is the complex amplitude of the velocity disturbance associated with the vortical 
excitation. The relative exit Mach number from the rotor can be related to the inlet Mach 
number of the stator through the relation M' +(X> = nooM-oo, since the nondimensional 
inlet mean velocity to the stator has magnitude one. The complex amplitude of the nth 
component of the disturbance entropy is given by 

^ — oo,n = (1 'T)'^-|-oo-^-oo w — oo,n (2.16) 

where the Mach numbers M' +00 and M_ 00 are based on the relative flow speeds, V ^ and 
TLoo, respectively. 

In summary, the wake flows are assumed to be parallel, aligned along the rotor exit flow 
angle, with constant total enthalpy and constant static pressure. At an axial location corre- 
sponding to the leading edge of the stator row, the wake velocity distribution is calculated 
using a shape function (cf. Eq. (2.7) or (2.8)), a centerline velocity defect at that location, 
which is found from Eq. (2.9), and a wake half- width at that location, which is evaluated 
from Eq. (2.10). By performing a Fourier decomposition on the component of the wake 
velocity normal to the inlet flow direction of the stator, values are found via Eq. (2.4). 
The magnitudes of entropic excitations corresponding to the values can then be deter- 
mined as well. The inlet vortical and entropic disturbances to the FEGV are characterized 
by interblade phase angles a n = na = tik^G, and frequencies u;„ = ruv = —nK v Vw, and are 
imposed at the inlet of the stator row, where an inviscid analysis numerically determines the 
resulting acoustic response. 

2.2 Unsteady Flow through the Stator Blade Row 

This section summarizes the method for calculating the unsteady pressure response of 
the FEGV to the imposed vortical and entropic excitations found from the wake analysis. 
The response is calculated using the LINFLO analysis (10, 11]. The entire analysis in this 
section will be performed in the stator frame of reference. 

The flow through the FEGV is considered to be a time-dependent, adiabatic, attached 
subsonic flow, with negligible body forces, of an inviscid non-heat-conducting, perfect gas. 
It is assumed that in the absence of unsteady excitation the mean flow is uniform at the inlet 
to the FEGV. The unsteady fluctuations in the flow arise from the entropic or vortical exci- 
tations at inlet, as well as the motions of control surfaces on the blades. All the excitations 
are assumed to be of small amplitude, periodic in time, and periodic in the 77 -direction. 
Control surface motions will be modeled as prescribed blade motions of the form 

'72- b ( x + tuG, t) = Re{r B (x) exp[7(u;f -f ma)\) x € B. (2.17) 

that are restricted to portions of the blade surfaces. Here "Rb measures the displacement 
of a point on a blade surface relative to its mean or steady-state position, x is a position 
vector, t is time, re is a complex displacement-amplitude vector, Re{ } denotes the real part 
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of { } and B denotes the reference (m = 0) blade surface. The incident disturbances are 
prescribed as functions of x and t which satisfy the field equations that govern the unsteady 
flow. Thus, small-amplitude entropic and vortical fluctuations relative to a uniform stream 
must be of the form 

5_oo(x,<) = Re{s_ 0O exp^K-oo • x + ut )]}... , £ < £- , (2.18) 

and 

C_<x>( x > t) = Re{<_oo exP^-oo -x + u>0]}--- , £ < • (2.19) 

Here s_ 0O and C-<x> are the (prescribed) complex amplitudes of the entropic and vortical 
fluctuations, s(x, t) and C(x, f), respectively, far upstream (f < £_) of the FEGV. 

The fluid motion is governed by a coupled set of nonlinear differential equations (i.e., 
the Euler equations) in continuous regions of the flow, a flow tangency condition at blade 
surfaces, and jump conditions at vortex-sheet unsteady wakes. In addition, information on 
the uniform flow conditions at inlet and exit and the entropic and vortical fluctuations at 
inlet must be specified. Since the unsteady excitations are assumed to be of small-amplitude, 
the time- dependent flow can be regarded as a small perturbation about an underlying mean 
or steady flow. This linearization can be performed by substituting relations analogous to 
that for the time-dependent fluid velocity, i.e., 

V(x,<) = V(x) + v(x, t) -F ... , (2.20) 

where V(x) is the local mean velocity, and v(x,f) is the first-order (in e) unsteady velocity 
into the full, time-dependent governing equations and collecting terms of like power in e 
and neglecting terms of 0(e 2 ) or higher. The first-order unsteady fluctuations that arise 
from the various independent modes of excitation are not coupled and hence, can be deter- 
mined separately. Indeed, it is sufficient to determine unsteady flow solutions for each single 
harmonic (in t and rj) component of a given disturbance. Solutions for arbitrary distur- 
bances and arbitrary combinations of various disturbances can then be obtained by Fourier 
superposition. 

The first-order or linearized unsteady flow properties caused by a periodic unsteady exci- 
tation at temporal frequency u> will be harmonic in time, e.g., v(x,f) = Re{v(x) exp{iu>t)}. 
In addition, the steady and, for an excitation at circumferential wave number = crG -1 , 
the first-order unsteady properties will satisfy the blade-to-blade periodicity conditions, e.g., 
V(x + mG) = V(x) and v(x + mG) = v(x) exp(im<r), respectively. Thus, solutions to time- 
independent nonlinear steady and linearized unsteady flow problems are required only over 
a single extended blade-passage region of the cascade. In addition, since analytic far-field 
solutions can be determined [14], the numerical solution domain can be restricted further to 
a single extended blade-passage region of finite extent in the axial-flow direction. 


The Steady Background Flow 


As a consequence of the assumption of uniform mean flow at inlet, the steady background 
flow will be isentropic and irrotational. Thus, V = V$, 


V-(pV$) = 0 


( 2 . 21 ) 
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and the fluid properties are related by 

=p^= = 1 - ^pMf^KV*) 2 - 1] , (2.22) 

where $, M , A, p and P are the steady velocity potential, Mach number, speed of sound 
propagation, density and pressure, respectively, and 7 is the specific heat ratio of the fluid. 
Numerical procedures for determining two-dimensional steady potential flows through cas- 
cades have been developed extensively, e.g., see [15] and [16], particularly for flows with 
subsonic relative inlet and exit Mach numbers (i.e., M Too < 1). The usual practice in such 
calculations is to solve the mass conservation equation (2.21), subject to a prescribed (i.e., 
VLoo and fl_oo) uniform inflow, flow tangency at blade surfaces, and a Kutta condition at 
blade trailing edges. 


The Linearized Unsteady Flow 

The system of field equations that govern the linearized unsteady perturbation of a 
potential mean flow can be cast in a very convenient form by introducing the velocity de- 
composition [17, 18] 

V = v, + V<f> = v R + V<f>* + V<t> , (2.23) 

where the rotational velocity, Vr, is taken to be divergence-free far upstream of the blade row, 
<f> , is a convected or pressure- less potential (i.e., D<f>*/Dt = 0) which satisfies the condition 
• n = — Vr • n at blade and wake mean positions, and the unsteady pressure depends 
only upon the potential, </>, through the relation p = —pD<f>/ Dt, where Dj Dt = iu> + V$ • V 
is a convective derivative based on the mean-flow velocity. 

The system of field equations that governs the linearized unsteady flow variables, s, Vr 
and <f> is determined by substituting the velocity decomposition (2.23) into the linearized 
Euler equations. If the underlying mean flow is isentropic and irrotational, the unsteady 
equations reduce to 


0 

II 

1 -KJ 

qIq 

1 

(2.24) 

-^(vr - sV$/2) + [(vr - sV$/2) • V]V$ = 0 

(2.25) 

and 


Wt (A ~ l! 5t ) “ f v ' 0*^) = f ' (**■) • 

(2.26) 

These equations are coupled only sequentially; hence, they can be solved in consecutively to 
determine the complex amplitudes of the entropy (s), rotational velocity (v/*) and velocity 
potential ( <f> ), respectively. 

Closed form solutions [17, 19] can be determined for the entropy and rotational velocity 
fluctuations in terms of the known conditions at inlet, i.e., 

s(x) = S-oo exp(iK_oo -X) , 

(2.27) 

and 

vr(x) = [V(X • A-oo) + s_ 00 V$/2]exp(tK_ 00 • X), 

(2.28) 
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where X = Aer + tfe/v, A and 'I' are the drift and stream functions, respectively, of the 
steady background flow, and A_oo = V/j,_ <» — 5_ 00 V_ 00 /2 [19]. Since C-oo = **-<» x v R,-oo 
and the rotational velocity is divergence-free far upstream of the blade row, the vectors /C-^ 
and vr-oo are orthogonal and Vr-oo = *( K - ~ x C-oo)/l #c -°o| 2 - 
A convected potential of the form [18, 19] 

4>. = [-iu-'A-oo ■ V-oo + F(#)] expiiK-n • X), (2.29) 


where 


f(*) = 


LO G COS — oo x A — oo) * e z 

2tt(1 — ia 0 u>) 


sin 


2tt[ , J , (x) - ^(x_)] 
Gcosfl-oo 


(2.30) 


will ensure that v, ■ n = (vr -f V<^„) • n = 0 at the mean blade and wake surfaces, and hence, 
that V<f> • n is finite at such surfaces. Thus, the unsteady flow variables s, Vr and <£* can be 
evaluated in terms of the drift and stream functions of the steady background flow. 

The unsteady velocity potential, <f>, is governed by the second-order partial differential 
equation, along with conditions at the blade and wake surfaces and at the inlet and exit 
boundaries. This function can be found numerically as a solution of the field equation (2.26) 
subject to the appropriate surface and far-field conditions. The flow tangency condition 


V<f> ■ n = [iu>r + (V$ • t)(t • V)r - (r • V)V$] • n (2.31) 


applies at the mean blade surfaces (i? m ), and the linearized unsteady pressure and normal 
velocity component are continuous across blade wakes ( VT m ), i.e., 

[D<f>/Dt\ = 0 and (V<£] • n = 0. (2.32) 


The velocity potential fluctuations in the far upstream and far downstream regions de- 
pend upon the prescribed excitation as well as on the acoustic and vortical response of the 
cascade. Analytic solutions for <j> can be determined [14], which satisfy the requirements that 
acoustic response disturbances either attenuate with increasing axial distance from the blade 
row or propagate carrying energy away from or parallel to the blade row and that vorticity 
must be convected downstream. These solutions can be matched to a near-field numerical 
solution and thereby serve to complete the boundary-value problem for the unsteady poten- 
tial. The vortical excitations determined by the wake analysis can thus be imposed at the 
inlet of a single extended blade passage of the FEGV, and the resultant acoustic response 
can be calculated. The amplitudes of the propagating acoustic response modes are then the 
targets of the active noise control analysis, which is presented in the next section. 
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3. Control Surface Analysis 


The analyses in the previous two sections described how vortical and entropic excitations, 
that represent the effects of the rotor wakes on the stator row, are determined and how those 
excitations are used in the LINFLO analysis to determine the subsequent acoustic response 
of the stator row. The propagating far-field acoustic response waves thus found make up 
the discrete tone noise at multiples of the blade passing frequency, and are the targets of 
the control analysis given in this section. Note that for the purposes of the control analysis, 
all of the acoustic response waves are considered to be independent modes. The means, for 
instance, that the upstream and downstream moving acoustic response waves at 2 BPF are 
considered different, even though they both have the same frequency and interblade phase 
angle. 

3.1 Anti-sound Generator Implementation 

Wake/blade-row interaction noise can be reduced by using discrete blade surface actuators 
to generate sound out of phase with the undesired acoustic responses. These actuators, or 
pistons, have been implemented into the LINFLO analysis as oscillating control surfaces 
on the suction and pressure surfaces of each blade. A general blade motion is represented 
by Eq. (2.17). Control surfaces, that oscillate with a prescribed translational displacement 
normal to the mean blade surface with a given frequency and interblade phase angle, are 
defined over restricted portions of the blade surfaces. For a set of K control surfaces or 
pistons, the complex amplitude of the displacement at the reference (m = 0) blade is given 

by 

K 

f b(x) = Er*[U(x- x Ue ) -U(x- x* ite )] (3.1) 

k=l 

where the index k ranges over the K pistons, r*. = r*ejv is the displacement of the kth piston 
with complex amplitude r*, is a unit vector normal to the blade surface, U is the unit 
step function, and x^je and Xfc >te are the locations of the leading and trailing edges of the 
fcth piston. Individual pistons are specified by choosing midpoint locations and lengths as a 
fraction of chord. The control analysis then calculates chordwise positions of the leading and 
trailing edges of the piston, and requires the surface segment within this range to perform a 
translational oscillation at the specified frequency with the specified complex amplitude. 

The analysis presented below shows how to determine the complex amplitudes r*, k = 

1 ,K of the pistons in order to best minimize the noise level, given the piston lengths and 
locations. The analysis assumes that there are K actuators on each blade, whose motions 
at each chordwise position are related to their counterparts on the reference blade through 
the interblade phase angle. An assessment of the effect of varying piston positions is given 
in the Results section. 

3.2 Control Algorithm 

The unsteady aerodynamic equations governing small unsteady perturbations of a uni- 
form mean flow can be solved to yield the wavenumbers and attenuation constants for all 
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acoustic response disturbances [14]. The form of the acoustic response, i.e., the wavenumbers, 
attenuation constants, and propagation directions of each response wave, is independent of 
the means of excitation. The linearity of the unsteady equations for small amplitude excita- 
tions implies that superposition is valid, so if both an external excitation and a set of control 
surface oscillations act simultaneously, the net complex amplitudes of the acoustic response 
waves will be equal to the sum of the responses to the individual excitations. Thus, the 
control excitations can in principle be “tuned” so that their corresponding acoustic response 
waves are of equal amplitude but 180 degrees out of phase with the undesired responses 
produced by vortical and entropic excitations. Such surface motions then act as sources of 
anti-sound which cancel the undesired noise. 

Assume that there are N acoustic response modes and K independent control surfaces. 
The total acoustic response due to both a specified excitation and a set of control surface 
motions will contain contributions attributable to each. The complex amplitudes of the net 
propagating acoustic response modes can therefore be expressed as a vector Pr such that 

Pr = Pn.exc + Pfl.ctri = Pr,«c + Ar . (3.2) 

The N x K matrix A contains as its elements the complex amplitudes of each propagating 
acoustic response wave resulting from a unit amplitude oscillation of each control surface, and 
the vector r contains the complex amplitudes of the control surface motions. If K = N, 
the matrix A in this equation is square and invertible. Consequently, a set of control surface 
amplitudes given by 

Texact = -A _1 pR, e xc ( 3 - 3 ) 

will yield p R = 0. In this case, the net amplitude of each propagating acoustic response 
wave will be zero. 

If, however, K < N, then the length of the total response vector ||pr|| = (Pr ■ Pn ) 1/2 
provides a natural measure of the net sound. In this expression, the superscript H denotes 
the Hermetian, or conjugate transpose, of the total response vector. In order to minimize 
||Pr|| in a least-squares sense, it can be shown [20] that it is necessary to chose the vector 
P«,ctri = Ar to be perpendicular to pr, i.e., 

Pfl.ctri ‘ Pfl = ( Ar) H ( Ar + PR, e xc) = 0 . (3-4) 

Expanding the Hermetian product allows this equation to be written as 

r H (A H Ar + A^Pr^xc) = 0 (3.5) 

whose nontrivial solution is given by 

riaq = — (A w A) -1 A w pR ie xc- ( 3 -6) 

This set of control surface motions will minimize the length of the acoustic response vector, 

IIphII- 

The least squares minimization treats each element of the propagating response vector 
as equally important. It may be true in certain situations that the reduction of some modes 
is more important than others. Consequently, a weighting matrix W can be introduced into 
Eq. (3.2), i.e., 

Wpft = Wp fl(IC + WAr. (3.7) 
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The matrix W is chosen to be square, diagonal, and composed of real elements. The values 
specified on the diagonal represent the relative importance assigned to the individual modes. 
The least-squares solution of this equation minimizes the length of a weighted net propagating 
response vector, ||Wp fi ||, and is given by 

r.l.q = -(A"W , WA)-’A"W I 'Wp J1 „ . (3.8) 

Note that if the matrix A is invertible, then Eq. (3.8) reduces to the solution given in 
Eq. (3.3), and that the least squares solution given in Eq. (3.6) represents the case where 
the weighting matrix is chosen to be the identity matrix. Note also that this procedure 
minimizes the length of the weighted net propagating response vector, ||Wp fl ||, which could 
result in a larger value for ||pfl|| itself than in the cases without control applied. 

The conclusion, therefore, is that if there are more propagating acoustic response modes 
than available independent control surfaces, a natural measure of the net sound produced 
is provided by a scaled length of the total acoustic response vector. The scaling is provided 
through the use of a specified weighting matrix W. The length of the scaled acoustic response 
vector is minimized by performing a least squares minimization procedure, whose result is 
given by Eq. (3.8), or by Eq. (3.6) if no weighting is performed. 


3.3 Circumferentially Averaged Sound Pressure Level 


The length of the total acoustic response vector ||pfl|| is directly related to the circum- 
ferentially averaged sound pressure level SPL. The sound pressure level, SPL, is defined to 
be 

SPL = 201o glo ^, (3.9) 

Pref 

where is the root mean square value of the unsteady pressure and p re f is a reference 
pressure, whose dimensional value is generally taken to be 2 x 10 -5 Pascals for airborne 
sound. 

The mean-square pressure, < p 2 >= p^, for a pressure disturbance at frequency u> is 
defined to be 

<P 2 >= hC P* d ( ut ) (3-10) 

The acoustic response at a single frequency is composed of a series of individual modes with 
different wavenumbers, i.e., 


p = Re{p exploit)} = Re ((! Pnj exp(iu><) j 


where the complex amplitude of the nth mode, p n , can be expressed as 

Pn = On exp[/3 n { + t(fc f ,„£ + *„,„!))) 


(3.11) 


(3.12) 


where a„,l3„, and n are the complex amplitude, attenuation constant, axial wavenum- 
ber, and circumferential wavenumber of the nth mode, respectively. Substituting Eq. (3.11) 
into Eq. (3.10) and performing the averaging results in 


<p 2 



+00 j 

£ (PmPn +PnP* m ) = 5 

n,m= — OO " 


+ OO 

E 

n,m=— 00 


(3.13) 
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where the superscript * denotes the complex conjugate. The terms in Eq. (3.13) depend 
upon circumferential position, by Eq. (3.12). Therefore, taking the square root of < p 2 > 
and substituting the result into Eq. (3.9) will result in a sound pressure level that is also 
dependent upon circumferential position. This dependence can be eliminated by averaging 
< p 2 > around the circumference, i.e., integrating from p = 0 to tj = NgG , where Ng is 
the number of blades in the stator row, dividing by the circumference, and substituting the 
resulting expression into Eq. (3.9). This result is defined to be the circumferentially averaged 
sound pressure level, SPL, which is given by 


/ +0 ° 1 \ 

SPL = 101og 10 I £ -p„p-J — 20 logjQ Pret 


(3.14) 


where the overbar indicates circumferential averaging. 

As above, consider the case where there are N propagating acoustic response modes. If 
the complex amplitudes of the propagating response modes are placed in a vector p«, then 
the argument to the logarithm in Eq. (3.14) is related to the magnitude of this vector via 

£ \pR,nPR.n = ipgPH = |||P«|P- (315) 

N prop ~ 

The logarithm is a monotonically increasing function, so the minimum of its argument cor- 
responds to the minimum of SPL. Consequently, since the least squares s olutio n for the 
control surface motions given in Eq. (3.6) minimizes ||p/i||, it also minimizes SPL. 

In summary, given N propagating acoustic response modes and K independent control 
surfaces, then if K = N the oscillation amplitudes of the pistons can be determined using 
Eq. (3.3) to give complete cancellation of the noise. If K < N , then the magnitude of the 
net propagating acoustic response vector provides a convenient measure of the net sound, 
and can be minimized by using the least squares minimization procedure given in Eq. (3.6) 
or the weighted least squares minimization procedure given in Eq. (3.8). The least squares 
minimization procedure also minimizes the circumferentially averaged sound pressure level, 

SPL. 
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4. Results 


The following results are based on an analysis of a representative fan stage composed of an 
upstream fan rotor and a downstream fan exit guide vane. The FEGV blades are constructed 
by superimposing the thickness distribution of a NACA 0008 airfoil on a circular arc camber 
line. The latter has a height at blade midchord of 0.08. The FEGV cascade has a gap to 
chord ratio G of 0.836, and a stagger angle 0 of 10 deg. The inlet and exit relative freestream 
Mach numbers are = 0.55 and M +00 = 0.45, respectively, and the inlet and exit flow 

angles are fi-oo = 30 deg and D +co = —2.39 deg respectively. The rotor blade passing 
frequency (BPF) as seen from the stator is lo = 4.106, and the interblade phase angle at 
BPF is a = —2.513. 

The chord length of the fan blades is approximately 2.3 times the chord of the FEGV 
blades. The gap to chord ratio G' of the fan cascade is 2.09. The wheel speed at the radial 
location chosen for the 2D cascade section is 1.366, nondimensionalized with respect to the 
inlet flow velocity of the response row. This set of parameters results in an exit flow angle 
for the rotor row of = —45 deg. 

4.1 Generation 

The wake model given in Section 2.1 was used to compute the velocity distribution down- 
stream of the rotor. A Fourier decomposition of this velocity field was then performed to 
determine the Fourier components of the unsteady excitation at the FEGV inlet. The mag- 
nitudes of the v+ components were evaluated at the leading edges of the stator blades. This 
location is two (rotor) blade chords downstream of the fan. Table 4.1 shows the chosen v+ 
values, along with the frequencies and interblade phase angles for the first three harmonics of 
the blade passing frequency. The analysis in Section 2.1 showed how the vortical disturbance 
amplitudes could be used to determine the amplitude of entropic disturbances as well, but 
the entropic disturbances at these amplitudes do not produce significant acoustic responses 
and are therefore not shown. 

An analysis of the equations governing the unsteady flow in the far-field of the stator 
reveals that at BPF, the acoustic response is “cut-off,” i.e., all of the acoustic waves present 
in the flow are of the decaying type. At 2 BPF, there are two propagating response waves, 


Vortical Excitation Data 


n 



< 

i 

-2.513 

4.106 

(.14807 x jO?*i0.) 

2 

-5.026 

8.212 

(.82745 x 10" 2 ,0.) 

3 

-7.539 

12.318 

(.41493 x 10- 2 ,0.) 


Table 4.1: Interblade phase angles, frequencies, and complex amplitudes for the first three 
harmonics of the vortical excitation calculated from the wake model. 
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Figure 4.1: Vortical excitation and acoustic response fields, (a) vortical field at BPF; acoustic 
fields at (b) BPF, (c) 2 BPF, and (d) 3 BPF. 


one in the upstream flow field and one in the downstream field, and at 3 BPF, there are four 
propagating response waves, with two in each flow field. 

Figure 4.1 shows contours of the real part of the unsteady vorticity field at BPF, and 
contours of the real part of the unsteady pressure at BPF, twice BPF, and three times 
BPF. The vortical fields at twice and three times BPF are similar to that at BPF, with 
slightly different shapes to the contours. These contour plots are the results of a numerical 
solution of the equations governing the steady and the first-order unsteady flows. In each 
case, vortical excitations have been imposed at the inlet to the cascade. The vorticity in 
Fig. 4.1a is convected through the blade row by the steady background flow. For the acoustic 
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response fields, the contour ranges were chosen to emphasize the propagating acoustic modes. 
Note that in the BPF acoustic contour plot shown in Fig. 4.1b, all the response modes are 
decaying, so the magnitude of the unsteady pressure drops quickly with increasing distance 
from the blade row. 

The pressure contours for the 2 BPF case in Fig. 4.1c show a single propagating wave 
upstream of the cascade and a single propagating wave downstream of the cascade. If all 
propagation angles are measured counterclockwise from the positive £ axis. The upstream 
acoustic wave is propagating at an angle of -170.8 deg and has a wavelength of 0.676 blade 
chords, while the downstream acoustic wave is propagating at an angle of -27.8 deg and has a 
wavelength of 1.98 blade chords. The complex amplitude of the upstream propagating wave 
is Pr,-oo = (—0.1309 x 10 -2 , 0.4937 x 10 -3 ), and the complex amplitude of the downstream 
propagating wave is Pr,+ <» = (0.1970 x 10 -2 , —0.2540 x 10 3 ). 

The 3 BPF acoustic field shows two propagating waves in each flow field. The upstream 
waves are propagating at angles of 172.8 deg and — 156.5 deg, and have wavelengths of 0.52 
and 0.42 blade chords, respectively. The downstream waves are propagating at angles of 
18.7 deg and -79.0 deg, and have wavelengths of 1.33 and 1.03 blade chords. The complex 
amplitudes of these propagating waves are 


Pr,- oo,i — (0.2065 x 

Pr,- oo ,2 = (-0.5222 x 

Pr, + oo,i = (0.1135 x 

PR, +00,2 = (0.8927 X 


10- 3 , 

0.4054 x 10" 3 ) 

10“ 3 , 

-0.3309 x 10" 3 ) 

10- 3 , 

0.2421 x 10" 3 ) 

10" 3 , 

0.2673 x 10" 3 ) 


(4.1) 


A linear superposition of two acoustic waves of equal magnitude and different propagation 
directions results in a series of ellipses on a contour plot of the unsteady pressure. The 
magnitudes of the upstream propagating waves at 3 BPF are relatively close to each other 
(a factor of approximately 1.36), and the angle between them is approximately 30 degrees. 
The acoustic field in Fig. 4. Id reflects this by showing somewhat elongated ellipses upstream 
of the cascade. The downstream flow field is dominated by the n — 2 mode, which is much 
larger in magnitude (by a factor of about 3.5) than the downstream n = 1 mode, and the 
contour plot shows this as well. 


4.2 Control 


Complete Cancellation 

The control procedures described in Section 3 have been implemented for the represen- 
tative fan stage. The acoustic responses at 2 BPF and 3 BPF contain propagating acoustic 
response waves with the complex amplitudes given above. To demonstrate the possibility 
of complete cancellation of the propagating acoustic modes, at each multiple of BPF the 
number of control surfaces was chosen to be equal to the known number of propagating 
acoustic response modes. For this demonstration, the lengths of all of the control surfaces 
was chosen to be 0.1 chord. The locations of the control surfaces were specified arbitrarily, 
and as long as non-overlapping control surfaces were selected, the resulting control matrix 
A in Eq. (3.3) could be inverted. The elements of A were determined by specifying a unit 
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Figure 4.2: Acoustic fields showing complete cancellation of propagating modes at (a) 2 BPF 
and (b) 3 BPF. 


amplitude displacement for each control surface as an excitation for the LINFLO analysis, 
and calculating the resulting complex amplitudes of each of the propagating response modes. 

Since two propagating acoustic response modes exist at 2 BPF, two control surfaces 
(on each blade) were selected. One piston was placed on the suction surface centered at 
0.2 chord; the other, on the pressure surface at 0.8 chord. The solution of Eq. (3.3) for 
the 2 BPF excitation gave the complex amplitudes of the control surface motions as n = 
(0.2542 x 10 -3 , -0.4153 x 10 -3 ) and r 2 = (-0.1224 x 10 -2 , 0.1080 x 10 -2 ). 

Four control surfaces were used for the 3 BPF case, two of which were placed on each 
surface of the blade. The control surfaces on the suction surface were centered at 0.2 and 0.7 
chord; those on the pressure surface, at 0.3 and 0.8 chord. The solution to Eq. (3.3) for the 
responses in the propagating modes to the vortical excitation and to unit amplitude piston 
oscillations yields 

n = (0.4286 x 10" 4 , 0.9962 x 10~ 4 ) 

r 2 = (0.9207 x 10“ 5 , 0.1902 x 10" 3 ) , , 

r 3 = (0.2665 x 10" 3 , 0.6581 x 10" 4 ) 1 ' } 

r 4 = (-0.4100 x lO" 4 , 0.1854 x 10~ 3 ) 

for the complex amplitudes of the control surface motions. The magnitudes of the acoustic 
response modes when both the vortical excitation and the control surface motions were 
present were all on the order of 10~ 7 . Contour plots of the acoustic fields with the control 
motions applied along with the vortical excitation are shown in Fig. 4.2. The response within 
the cascade in these figures is dramatically larger than that without the control motions, 
but the far field acoustic response has been virtually eliminated. No other structural or 
aerodynamic effects that may arise due to the piston motions have been evaluated. Should 
difficulties arise, however, the control analysis can be performed with the pistons at alternate 
locations. The effect of piston location on the control effectiveness is examined below. 
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Effect of Piston Location 


Since there are a finite number of grid nodes on the surfaces of the blades, the calculated 
leading and trailing edges of the piston did not fall directly upon grid points. The length 
of the piston as seen by the numerical analysis, therefore, was not exactly equivalent to the 
length specified, and a small change in the location of the midpoint of the piston for a given 
length yielded differing numbers of included gridpoints. This is an inevitable consequence 
of the non-uniform distribution of mesh points along the blade surface. The numerical 
length of the piston thus seen by the numerical code appeared to increase and decrease 
as it was moved along the blade, even though the specified length was not intended to 
change. The numerical representation of the piston thus implies that its specified length 
and location are not completely independent. This complicates parametric studies of the 
effects of changing lengths or locations of the control surfaces. Also to be noted is that the 
surface boundary condition contained in Eq. (2.31) contains the gradient of the displacement. 
The pistons have finite length, so their designated motions take the form of square pulses 
of translational oscillations of a given complex amplitude. The numerical implementation, 
however, performs a second-order central difference of the displacement on an irregular mesh, 
which has the effect of rounding the corners of the specified displacement. Neither of these 
effects, i.e., the weak coupling of piston length and location by the surface grid nodes, and 
the finite difference smoothing of the specified displacements, dramatically affect the control 
analysis. Rather, these factors are endemic to the numerical formulation, and make the 
results somewhat dependent on the density of the grid points on the surfaces of the blades. 
Each of the cases presented here was run on a global mesh containing 155 streamwise points 
and 40 circumferential points, which distributed 65 points on each of the suction and pressure 
surfaces of the blades. The number of surface grid points contained within the 0.1 length 
pistons ranged from four to seven, depending on location. 

To ascertain the effect of varying the control surface locations, a series of 50 tests was run 
at 3 BPF. In each case, a piston of length 0.1 was placed at varying locations on either the 
suction or pressure surface of each FEGV blade, and oscillated at unit amplitude. Twenty- 
five top and twenty-five bottom surface locations were examined. The magnitude of the 
acoustic response of each propagating mode is shown as a function of the piston location in 
Fig. 4.3. Fig. 4.3a contains information from moving the piston along the suction surface, 
while Fig. 4.3b contains information from moving the piston along the pressure surface. 
There are four propagating modes at this frequency and interblade phase angle, so four lines 
appear on each figure. 

An examination of Fig. 4.3 suggests several conclusions. First, it appears that there 
are no obviously preferred locations for the pistons. Each mode has a completely different 
reaction to variations in piston location, and while the curves do show maxima, the maximum 
response value of each curve occurs at a different piston location. Second, the figure shows a 
certain jaggedness in some of the response magnitude curves. This effect is due to the weak 
coupling between length and location described above. Third, since each propagating mode 
has a unique dependence upon piston location, any conclusions made on piston location 
would be case dependent. While a complex optimization could in principle be performed 
to find a set of piston locations that would minimize the necessary control amplitudes, it is 
unlikely that the resulting set of locations could be generalized to arbitrary cases. Finally, 
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Figure 4.3: Effect of piston placement on mode magnitudes, (a) Suction surface piston, (b) 
Pressure surface piston. In each case, the solid line is for mode 1 upstream, the short dashes are 
for mode 2 upstream, the mixed long and short dashes is for mode 1 downstream, and the long 
dashes are for mode 2 downstream. 


the acoustic response modes generally have only a limited dependence on actuator location. 
The greatest dependence observed in this study occurred for the n = 1 upstream mode, 
which varies in magnitude by only a factor of three when the piston is moved along the 
suction surface. It would seem, therefore, that no “rule of thumb” conclusions about piston 
placement can be drawn, and that while an optimization of locations could be performed, 
the resulting improvement in control displacements is not likely to be dramatic. 
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Least Squares Minimization 


To illustrate the effects of the least squares minimization solution given in Eq. (3.6) and 
the weighted least squares solution given in Eq. (3.8), consider the same 3 BPF case as above, 
which contains four propagating acoustic response waves. Assume now that there are only 
two pistons available for control. The two pistons selected for demonstration of the least 
squares concept are the ones used to achieve complete cancellation for the 2 BPF excitation. 
Recall that in that case, the pistons each had length 0.1, and were centered at x = 0.2 on 
the suction surface and x — 0.8 on the pressure surface, respectively. 

The least squares minimization is applied by solving Eq. (3.6), which yields the complex 
amplitudes n = (-0.4830 x 10“ 5 , -0.4022 x 10" 4 ) and r 2 = (-0.7891 x 10‘ 4 , -0.8567 x 10" 5 ). 
The acoustic contours for this case are shown in Fig. 4.4a. With the control applied, the 
changes in the magnitudes of each propagating response wave are shown in Table 4.2. The 
magnitude of the net propagating response vector \\pr\ \ has been reduced by approximately 
28.4% overall, which corresponds to a decrease in the circumferentially averaged SPL of 
about 3 dB. 

To demonstrate the weighted least squares minimization procedure, a weighting matrix 
is chosen to require complete cancellation of both downstream modes. The required control 
displacements for this case are then r q = (—0.5216 x 10 -3 , 0.3196 x 10~ 3 ) and r 2 = (—0.2631 x 
10 -3 , 0.2343 x 10“ 3 ). Acoustic contours for this case are shown in Fig. 4.4b. Note that the 
downstream flow field is completely quiet, while the upstream response has been increased 
significantly. The results are summarized in Table 4.3. 

These results are dramatic, in that, although the downstream acoustic response waves 
have been virtually eliminated, the magnitudes of the upstream response waves have in- 
creased by almost two orders of magnitude. An examination of Fig. 4.3 reveals the reason 
for the strong upstream response. In the figure, it is clear that the downstream acoustic 
response waves have uniformly lower magnitudes than the upstream response waves, when 
excited by a unit amplitude control excitation. Consequently, a relatively large control 
displacement is required to achieve a large enough response to cancel the undesired down- 
stream modes. This large control surface displacement is then multiplied by the far greater 
sensitivity of the upstream modes to yield a strong upstream response. 


Leasts Squares Minimization 
3 BPF with 2 Pistons 

Magnitudes of Each Propagating Acoustic Response Mode, Hp^nll 


n 

No control 

Control 

Percent Change 

i 

0.4550 x 10“ 3 

0.2738 x 10~ 3 

-39.8% 

2 

0.6182 x 10~ 3 

0.3931 x 10“ 3 

-36.4% 

3 

0.2674 x lO" 3 

0.2872 x 10" 3 

+7.40% 

4 

0.9319 x 10~ 3 

0.6876 x 10" 3 

-26.2% 


Table 4.2: Least squares minimization of the net acoustic response, with four propagating 
acoustic response modes and only two control surfaces. 
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a) Least Squares Minimization 


bi Downstream cancellation onl 
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c) Upstream cancellation onl 
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Figure 4.4: Weighted least squares minimizations, (a) Least squares solution with two pistons 
and four modes, (b) Downstream cancellation only, (c) Upstream cancellation only. 

This suggests that if the weighting matrix is chosen so that only the upstream modes 
are cancelled rather than the downstream modes, the overall result will not be as severe. 
Application of the control analysis reveals that this is indeed the case, and acoustic contours 
for this case are shown in Fig. 4.4c. The required control displacements are rq = (-0.1293 x 
10~ 4 , -0.3205 x 10 -4 ) and r 2 = (-0.3094 x 10 -4 , 0.3616 x 10 -4 ), and the changes in the 
response modes are shown in Table 4.4. The upstream response modes respond more readily 
to the control, and this is reflected in the fact that even though the downstream modes were 
ignored in the weighted analysis in this case, one of them actually decreased in magnitude 
by almost 25%. 

The conclusion to be drawn, therefore, is that in each case where control is to be applied 
and there are not enough pistons available for complete cancellation, it is important to 
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Weighted Leasts Squares Minimization 

3 BPF with 2 Pistons 
Downstream Cancellation Only 

Magnitudes of Each Propagating Acoustic Response Mode, ||pR,n|| 


n 

No control 

Control 

Percent Change 

i 

0.4550 x 10~ 3 

0.1404 x 10" 1 

+2986% 

2 

0.6182 x 10" 3 

0.8098 x lO" 2 

+1210% 

3 

0.2674 x 10~ 3 

0.2066 x 10" 6 

-99.92% 

4 

0.9319 x 10~ 3 

0.3232 x lO” 6 

-99.97% 


Table 4.3: Weighted least squares minimization cancelling the downstream acoustic response 
modes. 


Weighted Leasts Squares Minimization 

3 BPF with 2 Pistons 
Upstream Cancellation Only 

Magnitudes of Each Propagating Acoustic Response Mode, ||pH,n|| 


n 

No control 

Control 

Percent Change 

1 

0.4550 x 10~ 3 

0.1165 x 10~ 6 

-99.95% 

2 

0.6182 x 10- 3 

0.2156 x 10~ 6 

-99.97% 

3 

0.2674 x 10“ 3 

0.2001 x 10" 3 

-24.98% 

4 

0.9319 x 10" 3 

0.9989 x 10- 3 

+7.190% 


Table 4.4: Weighted least squares minimization cancelling the upstream acoustic response 
modes. 

examine the behavior of each propagating acoustic response mode to the control. 
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5. Conclusions 


This report presents a combined analytical/computational investigation of sound gener- 
ation in cascades, and the subsequent active control of the acoustic responses through the 
use of blade-surface actuators. To ascertain the proper dimensional magnitudes of represen- 
tative acoustic response, an existing fan wake model was analyzed to determine appropriate 
unsteady vortical and entropic excitations to a downstream FEGV. The magnitudes of the 
excitations depend upon the geometry of the fan stage, the wheel speed and exit mean flow 
of the rotor, and the shape, centerline velocity defect, and half-width of the rotor wakes. 
Examples of the generation of acoustic response waves were presented by applying vorti- 
cal excitations at the inlet of the FEGV, and solving the governing steady and first-order 
unsteady flow equations numerically using the LINFLO analysis. 

At a given frequency and interblade phase angle the form of the acoustic response of 
the FEGV, i.e., the wavenumbers, attenuation constants, and propagation directions of the 
far-held acoustic response waves, is independent of the means of excitation. Therefore, the 
acoustic responses due to a set of control surface oscillations on the FEGV blades are of the 
same form as those due to an external excitation. The complex amplitudes of the control 
surface motions can be related to the complex amplitudes of the undesired propagating 
acoustic response waves through a set of simultaneous linear equations. If the number of 
available control surfaces is equal to the number of propagating acoustic response waves, this 
system of equations can be solved exactly, resulting in a set of control displacements that 
completely cancel the noise. Examples of complete cancellation at twice and three times 
BPF were shown. 

If the number of acoustic response waves exceeds the number of control surfaces, then 
the magnitude of a net propagating acoustic response vector provides a convenient measure 
of the effectiveness of the control. A least squares minimization procedure minimizes this 
quantity, which is directly related to the circumferentially averaged sound pressure level. 
A weighting matrix has also been incorporated into the control analysis, to allow selective 
cancellation of the acoustic response waves. An example of the least squares minimization 
procedure was presented, as well as cases where the noise was completely eliminated in either 
the upstream or the downstream flowfield of the FEGV . 

The magnitude of each propagating acoustic response wave has a unique dependence 
upon piston location. Consequently, no clear “rule of thumb” exists for piston placement, 
which is therefore case dependent. A degree of relative insensitivity of some of the modes 
to piston location was also observed in the cases studied, which implies that a complex 
optimization procedure to determine the piston locations that lead to minimum required 
control displacents is not likely to result in a dramatic improvement. 
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List of Symbols 


All physical parameters listed below are dimensionless. Lengths have been scaled with 
respect to blade chord, time with respect to the ratio of blade chord to upstream freestream 
flow speed, density with respect to its upstream freestream value, velocity with respect to 
the upstream freestream flow speed, and entropy with respect to the fluid specific heat at 
constant pressure. The number(s) in parentheses at the end of each symbol description 
indicates an equation in which the symbol appears. 


Roman 


•A— oo 

Constant vector, (2.28). 

A 

Speed of sound propagation in steady background flow. 

A 

Coefficient matrix in control equations. 

a 0 

Constant, (2.30). 

B 

Moving blade surface. 

B 

Blade mean position. 

Cd 

Section drag coefficient, (2.9). 

e 

Unit vector. 

F 

Complex function of mean-flow stream function, (2.30). 

f(N') 

Wake profile shape function, (2.7). 

G 

Cascade gap vector (= Ge n ), Fig. 2.1. 

H t 

Total enthalpy, (2.13). 

H t 

Total enthalpy in steady background flow. 

H 

Static enthalpy, (2.13) 

H 

Static enthalpy in steady background flow. 

i 

Imaginary unit. 

K 

Number of control surfaces per blade. 
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k{,k v Axial and circumferential wave numbers. 

M Mach number in steady background flow, 

m Blade number index 

P Pressure. 

P Pressure in steady background flow. 

p First-order unsteady pressure. 

p Complex amplitude of first-harmonic unsteady pressure. 

p R Vector containing complex amplitudes of the propagating 

acoustic response waves. 

Prm» Root mean square pressure, (3.9) 

H Surface (blade or wake) displacement vector, (2.17). 

r Complex amplitude of surface (blade or wake) displacement 

vector, (2.17); vector of complex amplitudes of actuator displacements, (3.1). 

r* Complex amplitude of displacement of fcth actuator. 

3 First-order unsteady entropy, (2.14). 

s Complex amplitude of first-harmonic unsteady entropy, (2.16). 

SPL Sound pressure level, (3.9). 

t Time. 

( T',N ') Coordinates tangential and normal to the exit 

flow direction of the rotor. 

u+,u + Components of perturbation velocity in Fig. 2.2 tangential 

and normal to the inlet mean flow velocity of the stator. 

V Velocity. 

V Velocity in steady background flow, 

v First-order unsteady velocity. 
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Vw 


V M 


V 


Vr 


V* 

w 

W 

W 

X 

x 


x,y 


Greek 


/? 

7 

A 

6 

c 

c 

0 

K 


Wheel velocity of rotor, = V^e^, Fig. 2.1 
Wake velocity defect, (2.6) 

Magnitude of wake edge velocity, (2.11) 

Minimum wake velocity, (2.6) 

Complex amplitude of first-harmonic unsteady velocity. 

Complex amplitude of first-harmonic rotational unsteady velocity. 
Complex amplitude of first-harmonic source-term unsteady velocity. 
Moving wake surface. 

Wake mean position. 

Weighting matrix for control equations, (3.7) 

Lagrangian coordinate vector. 

Position vector. 

Cartesian coordinates along and normal to mean position of blade chord, 
Fig. 2.1. 


Acoustic attenuation constant. 

Fluid specific heat ratio. 

Drift function. 

Wake half-width, (2.10) 

First-order unsteady vorticity. 

Complex amplitude of first-harmonic unsteady vorticity. 
Cascade stagger angle, Fig. 2.1. 

Wave number vector. 
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t/ Cascade axial and circumferential Cartesian coordinates, Fig. 2.1. 

p Density. 

p Density in steady background flow. 

p First-order unsteady density, 

er Interblade phase angle, 

r Unit tangent vector, (2.31). 

$ Velocity potential for inviscid steady background flow. 

^ Velocity potential for first-order inviscid unsteady flow. 

<f> Complex amplitude of first-harmonic unsteady velocity potential. 

<f> m Complex amplitude of first- harmonic convected potential. 

Stream function for inviscid steady background flow. 

D Steady flow angle, Fig, 2.1. 

ur Temporal frequency. 

Mathematical 


D/Dt 

D/Dt 

Im{} 

Re{} 

V 

[] 

< > 

INI 


Convective derivative operator. 

Convective derivative operator based on mean-flow velocity. 
Imaginary part of { }. 

Real part of { }. 

Gradient operator. 

Change in a flow quantity across a surface of discontinuity, (2.32). 
Time average, (3.10) 

Magnitude of a vector. 
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Subscripts 


B 

exc, R, ctrl 

JV, T, z, Tj 
and £ 

~ 5 + 

=foo 

Superscripts 

f 


H 

* 


Reference blade surface. 

Excitation, response, control. 

Indicates direction: N — normal to inlet freestream direction, 

T — parallel to inlet freestream direction, z — out from the 
page, tj — cascade “circumferential” direction, and 
£ — cascade axial direction. 

Axial locations at finite distance upstream and downstream from blade 
row; point of intersection (£_,77_) of axial line £ = £_ and 
reference blade stagnation streamline. 

Far upstream/downstream freestream value of a steady flow variable, 
Fig. 2.1; far upstream/downstream value of an unsteady flow variable. 


Quantity associated with rotor. 

Steady background flow quantity; 
Circumferential average (3.14) 

Time-dependent flow variable. 

Hermitian (conjugate-transpose) of a matrix. 

Complex conjugate, (3.13). 
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